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Abstract 



r^ . The incorporation of realistic trajectories of two-particle scattering into 

existing BUU or QMD codes meets considerable difficulties. We pro- 
pose a method of replacing the scattering event by a tractable non-local 
scenario reproducing the correct asymptotics. The first principle deriva- 
^ \ tion justifies the use of corresponding nonlocal kinetic equations. The 

^-v I nonlocal shifts are necessary if the system should approach an equation 

f^ ■ of state with second quantum virial coefficient. The modifications of 

^^ . QMD and BUU codes are discussed and comparison are presented with 

Q>^ \ resent data. We find a pre-equilibrium production of particles leading 

On ■ to higher energetic distributions of protons. 



1 Introduction 

Recent dynamical simulations of heavy ion reactions at low and mid energies 
are based either on BUU or QMD simulations. These equations simulate the 
Boltzmann equation (BE) with the extension of Pauli-blocking and meanfield 
drift 

dfi deidfi deidfi ^^ f dpdg 1 da < _ _ \ 

dt dk dr dr dk ~^J {2TrnYfi^dn^^^ ^^ ^' ^^' 

fMl - /i)(l - /2) - (1 - /3)(1 - /4)/l/2] . (1) 

The arguments of distributions / and energies e are shortened as /i = fa{k, r, t), 
/2 = fb{p,r,t), /s = faik-q,r,t), and f 4 = fbip + q, r, t), withmomenta k,p, q, 
coordinate r, time t, spin and isospin a, b and reduced mass /x. 

While the motion due to nonlinear meanfield is governed by the quasiclas- 
sical Hamilton equations for chosen test particles or wave packets, the collision 
is performed randomly if the particles meet inside the area of cross section or 
at the closest approach. By random choice of transfer momenta and averaging 



over many runs the scattering probability given by the cross section in the 
Boltzmann equation is reproduced. Medium effects are discussed in terms of 
density and temperature dependent depletion of cross section. 

However, there are modifications of the Boltzmann equation which have 
a more fundamental meaning. It has been noticed that the starting Boltz- 
mann equation has an inherent contradiction. While the collision probability 
describes a finite diameter the collision takes place at a time instant and local 
in space. This deficiency has to be cured. Simple cases like hard sphere gases 
have been discussed already be Enskog [jl| and in nuclear matter by Malfliet 
[|] and Halpert 0. Resonance scattering has been discussed by Pratt and 
Danielewicz |^ . The aim of these corrections are to describe a realistic scat- 
tering event within kinetic theory in a way that the correct asymptotics are 
reproduced. Only by this way the virial corrections, e.g. the two-particle cor- 
relation energy in balance equations can be expected to be approached during 
time evolution. The original Boltzmann equation and consequently the stan- 
dard BUU and QMD simulations lead to no binary correlation energy. As it 
was shown this binary correlation can be reproduced within a non-local 
collision integral. 

The need for nonlocal corrections can be stimulated by discussing the scat- 
tering of two particles as superpositions of wave packets and similar [0. 
The solution of the free Schroedinger equation can be written for the relative 
motion with cm. momentum k as 

0(p,fc,t) = F(p,fc)e-^^-* (2) 

and the asymptotic wave after scattering can be decomposed into \E' = + (p'^^ 
with the scattered wave for large distance x from scattering center 



^"<=^'*''' = /(2If^'^-^> ^'''T''"" °"""'" 



(3) 



with the scattering amplitude f{p, cos(px)). We proceed now and expand the 
scattering amplitude around the center of mass momentum k 

f{p,cos{px)) = |/(p,cos(px))|e''^(P'™^(''")) 

= f{k, cos{kx)) {l + {p- k)Vp\f{p, cos{px))\j,=k) e*{p-'=)Vp<5(p,cos(px)),.,_ 

(4) 

The derivative of the phase 5 leads now to the definition of the effective space 
shifts A and the time shift A|| 

o - cos(xA;)t|7 k 

2A = Vp6 = dk6 + ^ ^9co.(.fc)5 = -A|| + A^ (5) 

\k\ m 



where we denoted the shifts corresponding to the direction of k as || and _L. 
Expanding in (^) also Ep = p^/m = k'^/m + (p — k)2k/m and collecting all 
terms we obtain 

where the wave packet velocity is u = 2k/m. We observe three effects of 
scattering on the asymptotics: (i) a genuine time delay Ay, (ii) an effective 
displacement of the two colliding particle of A with respect to the center 
of mass and (iii) a modification of scattering probability JF = JF(1 + (p — 
fc)Vp|/(p,cos(px))|p=fc). 

In this contribution we will put these ideas of nonlocalities on the firm 
ground using the quantum kinetic equation with nonlocal scattering integrals 
which was derived from quantum statistics |]5|, |^ to show how the effect of 
nonlocalities play a role in simulations of heavy ion reactions and compare 
them with experiment. 

2 Nonlocal kinetic equation 

In p, 1^, 1^ we have given a systematic derivation of nonlocal kinetic equation 
(NKE) which includes hard sphere like corrections, resonance like corrections 
as well as rotational angle corrections, which turned out to be important for 
angular momentum conservation. This kinetic equation completes all balance 
equatons on the level of quantum second virial coefficients. The NKE-equation 
with the collected gradient terms reads [A^ = \{^2 + A3 + A4)] 

dfi ^ dei_df^ _ de^dfi_ ^y f dpdq ^ + ^^ - £3 - £4 + 2Ae) 



dt dk dr dr dk '^ J (27r) 
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(7) 



p+q 




Figure 1: The equivalent 
scattering event of two col- 
liding particles. 



with the nonlocal shifts given as derivatives of the total phase shift of the 
scattering T-niatrix < 
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After derivatives, A's are evaluated at the energy shell Q ^ e^ + £4. Unlike 
in (P, the subscripts denote shifted arguments: /i = fa{k,r,t), /2 = fb{p,r — 
A2,t), h = fa{k-q-AK,r-As,t-At), and U = fbip+q-AK,r-A^,t-At). 

The A's are effective shifts and they represent mean values of various nonlo- 
calities of the scattering integral as demonstrated in figure ^ These shifts enter 
the scattering integral in the form known from the theory of gases [jl], ^, |T0| , 
however, the set of shifts is larger due to the medium effects on the binary 
collision that are dominated by the Pauli blocking of the internal states of the 
collision. 

Now we want to use these derived nonlocalities to mimick the real scat- 
tering event in simulations. To this end we use a classical analogy of these 
nonlocalities and the scattering event as demonstrated by picture |l|. We ob- 
tain the following scenario: Two particles approach until they reach a distance 
A2. Then they form a molecule traveling over a distance Aj with a time Aj. 
During this propagation the molecule rotate about A^. Now the kinetic equa- 
tion (0) is numerically tractable by recent Monte Carlo or BUU codes [|l^ as 
demonstrated below. 
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Figure 2: A nonlocal bi- 
nary collision (solid line) to- 
gether with the scenario of 
sudden jump at the closest 
approach. 



3 Instant approximation 

The selfconsistent evaluation of all A's for all collisions would be too demand- 
ing. We employ two kinds of additional approximations. First, following ap- 
proximations used within the BUU equation, we neglect the medium effect on 
binary collision, i.e., use the well known free-space T- matrix. Second, we rear- 
range the scattering integral into an instant but non-local form. This instant 
form parallels hard-sphere-like collisions what allow us to employ computa- 
tional methods developed within the theory of gases |jl2| similarly as it has 
been done in [|lT . 



In the instant approximation we let particles to make a sudden jump at 
time t from r^ and r;, to effective final coordinates Va and ff,. These effective 
coordinates and momenta k and K are selected so that at time t + At particles 
arrive at the correct coordinates, r^ and r^, with the correct momenta, k' and 
K'. Accordingly, in the asymptotic region, after t + At, there is no distinction 
between the non-instant and instant pictures, which is shown as solid line in 
figure ^ This asymptotic condition is naturally met if one extrapolates the 
out-going trajectories from known coordinates and momenta at t + A^ back to 
the time t. Doing so one finds that the effective coordinates read 
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(9) 
(10) 



When incorporating the displacements into the QMD simulation code, we 
have to face the fact that two particles are selected for a collision if they meet 



at the point of closest approach. This distance is different from the distance A2 
required from the equivalent scattering scenario presented in figure || as solid 
line. We consider now the time required to travel from A2 to the distance 
of closest approach At = ^kA2 in analogy to |jl3|. Within this scenario we 
are allowed to jump at the point of closest approach to the final asymptotics 
(K) and ([TO|)with the additional distance the particle travel during A^. The 
effective final coordinates thus have to be evaluated as 



^^L±^^A, (11) 



with the effective displacement 



A = Ias - A3 + ^^(A, - At). (12) 

2 m 

Since the center of mass does not jump in the collision, the final displacement 
can be also written in an alternative way, A = |A2 + A4 — A2 — ^^{At — Aj). 
The non-local corrections are thus performed as follows. When the collision is 
selected, we evaluate A from (|12D and (||), redisplay particles into f^ and f^ 
and continue with the simulation. 

Let us note that from (H) and ([T^) the already presented shift (^ follows 



as it should. 

At this point it is possible to establish a connection of the present theory 
to the hard-sphere-like corrections used by Malfliet and Kortemeyer, Baffin 
and Bauer |[Tl|. For hard spheres of the diameter d, the phase shift has a 
classical limit (p = ^ ~ kM which gives A3 = and A2 = A4 = Ad. The 
displacement thus has the same amplitude d for all binary collisions and points 
in the direction of the transferred momentum, as it is known from the Enskog 
equation |[I|]. 

4 Results for Simulation 

In order to investigate the effect of non-local shifts on realistic simulations of 
a heavy ion reaction, we have evaluated A^ from the two-particle scattering 
T-matrix T^ in the Bethe-Goldstone approximation [iMl B using the separable 



Paris potential |T^ . The comparison of the shifts calculated for different po- 



tentials concerning partial wave coupling up to B-waves can be found in [|T6 
We have incorporated these shifts into a QMB code for the central collision of 
i^sXe^^i^Sn at 50 MeV/A. 

Figure |]a shows the exclusive proton spectra subtracting the protons bound 
in clusters. This procedure is performed within a spanning tree model which is 
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Figure 3: The particle spectra for central collision of ^^^Xe— i>^^^Sn at 50 MeV/A 
with and without non-local corrections. The data are extracted from recent 
INDRA experiments [1^. The non-local corrections bring the spectrum of the 
protons towards the experimental values leaving the clusters almost unchanged. 
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Figure 4: The number of 
collisions per time with and 
without non-local collisions 
within a BUU simulation of 
the same reaction as in fig- 
ure (i. 



known to describe a production of light charged cluster in a reasonable agree- 
ment with the experimental data, Figs. ^-f. Within the local approximation, 
however, the remaining distribution of high-energy protons is too low to meet 
the experimental values. As one can see, the inclusion of non-local collisions 
corrects this shortage of the QMD simulation. As demonstrated in Fig. ^, 
productions of light clusters are rather insensitive to the non-local corrections. 
This also shows that the improvement of the proton production is not on cost 
of worse results in other spectra. 

A microscopic mechanism leading to the increase in the high-energy part of 
the particle spectrum can be traced down to an enhancement of the number of 
collisions at the pre-equilibrium stage of the heavy ion reaction demonstrated 
in Fig. ^ for the BUU simulation of the same reaction. This enhancement 
gives rise to an immediate proton production which itself translates into a 
high energetic spectra. In other words, the strong production of the high- 
energy protons follows from the pre-equilibrium emission of particles. The 
BUU simulation also shows that non-local corrections are important namely 
in the early stage of reaction well before most of light clusters are formed. It 
explains why the production of protons is affected while the formation of light 
clusters is nearly untouched by the non-local corrections. 



5 Conclusion 

In summary, as documented by the improvement of the high-energy proton 
production, the non-local treatment of the binary collisions brings a desirable 
contribution to the dynamics of heavy ion reactions. According to an experi- 



ence from the theory of gases, one can also expect a vital role of non-localities 
in the search for the equation of state of the nuclear matter. It is encouraging 
that the non-local corrections are easily incorporated into the BUU and QMD 
simulation codes and do not increase computational time. Corresponding pro- 
grams can be obtained from authors. 
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